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In this talk 0, we present a recently suggested way on how to analytically incorporate massive threshold effects into observables 
calculated in massless QCD. No matching is required since the renormalization scale is in this approach connected to the physical 
momentum transfer between static quarks in a color singlet state. We discuss massive fermionic corrections to the heavy quark 
potential through two loops. The calculation uses a mixed approach of analytical, computer-algebraic and numerical tools including 
Monte Carlo integration of finite terms. Strong consistency checks are performed by ensuring the proper cancellation of all non-local 
divergences by the appropriate counterterms and by comparing with the massless limit. The size of the effect for the (gauge invariant) 
fermionic part of ay(q 2 ,m 2 ) relative to the massless case at the charm and bottom flavor thresholds is found to be of order 33%. 
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1 Introduction 

1.1 Analytic Thresholds at one Loop 

In the MS and MS renormalization schemes, the run- 
ning of the QCD coupling a s , by construction, does not 
"know" about the masses, m q , of quarks. The (3 function 
Elu describes the evolution of the strong coupling "con- 
stant" in the asymptotic regime, i.e. for values of the 
renormalization scale y, 3> m q . Near the quark flavor 
thresholds one has to turn to effective descriptions which 
match threories with n massless quarks onto a theory 
with n—1 massless and one massive flavor at the "heavy" 
quark threshold. In this way the dependence on the di- 
mensional regularization mass parameter \x is reduced to 
next to leading order effects by giving .un-the analyticity 
of the coupling at the flavor thresholdB'QB. 

While this procedure of matching conditions and ef- 
fective descriptions is certainly workable, from a theo- 
retical standpoint it would be advantageous to have a 
physical coupling constant definition which is analytic at 
thresholds. In addition, as a physical observable, the to- 
tal derivative with respect to the renormalization scale 
\x vanishes. Such a system is given by identifying the 
ground state energy of the vacuum expectation value 
of the Wilson loop as the potential V between a static 
quark-antiquark pair in a color singlet state bQqBEJO: 
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in momentum space with q 2 
factor Cp is the value of the Casimir operator T a T a in 
the fundamental representation of the external sources 
and factors out to all orders in perturbation theory. As 
one is free to choose the representation of the external 
particles, we obtain the static gluino potential by adopt- 
ing the adjoint-representation. 

Recently 113, the effect of the massive fermionic one 
loop contributions to the heavy quark potential were in- 
corporated into a continuous and smooth function rif (p q ) 
given to lowest order by 
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-q 2 . The mass dependence 



of the physical ^-scheme can now be transferred to the 
MS scheme by using the commensurate scale relation 
E£l between the tjsip schemes. Employing the multi scale 



approach of Ref . liH 
through two loops E 



cives the following scale fixed relation 



V{r,m 2 q ) = - lim hog(0\Tr {P exp ( j> dx„A£T a j }|0) 

(1) 

where r denotes the relative distance between the heavy 
quarks, m q the mass of "light" quarks contributing 
through loop effects and T a the generators of the gauge 
group. It is then convenient to define the effective charge 
ay{(l 2 ,m 2 ) as 
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Q* = 2.3Q , Q** = 6.539Q (5) 

whereas Q*** to this order is not constrained. A first 
approximation is obtained, however, by setting Q*** — 
Q** Ej. Note that the scale Q at one loop is a factor 0.4 
smaller than the physical scale Q*. 



1 




la r ~a rM ]/3 r 

0.004 



Q [GeV] 




Figure 1: The sum over massive quark flavors for the effective 
analytic function h°—— given in Eq. ^. 



1.2 Analytic a 



MS 



One is now free to adopt the commensurate scale relation 
of Eq. |^ as a definition of the extended scheme 6ej^t3. 
At one loop we have therefore 



a Trs(Q) = a v{Q* 
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for all scales Q. Eq. |TJ not only provides an analytic ex- 
tension of MS like renormalized schemes, but it also ties 
down the renormalization scale fj, to the physical scale 
with massive quarks, entering into the vacuum polariza- 
tion contributions to ay . There is thus no scale ambigu- 
ity in perturbative expansions in ay or o.jjg- To lowest 
order we obtain in addition 
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A very good (~ 1%) approximation is given by the fol- 
lowing simple result for the one loop effective function of 
flavors: 
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Fig. [TJ contains the analytic function nf summed over 
all quark flavors in comparison with the conventional 
step-function approach, which treats quarks as infinitely 
heavy below threshold and massless above. 



Figure 2: The relative difference between the calculation of a^Y^(s) 

in the analytic extension of the MS scheme and the standard treat- 
ment of masses in the MS scheme. The discontinuities are due to 
the mismatch between the s/ m y | jir| d m 2 / s expansions of the ex- 
plicit QCD corrections in Refs. 



1.3 Applications 

It is of course possibJe to treat mass effects exactly within 
the MS schemecJ'Ej by explicitly calculating higher twist 
QCD corrections. In order to have a meaningful com- 
parison between the different approaches of incorporat- 
ing mass effects, we choose to compare the above treat- 
ment to the explicit corrections to an observable, here 
the quark part of the non-singlet hadronic width of the 
Z-boson, . Writing the QCD corrections in terms 

of an effective charge we have: 
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where the effective charge a^ g (s) contains all perturba- 
tive QCD corrections. At the one loop level we are left 
with a very simple expression 
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This simple expression reflects the fact that the effects 
of quarks in the perturbative coefficients, both massless 
and massive, should be absorbed into the running of the 
coupling. The BLM-scale iA is Q* = 0.7076^5 for this 
observable E3. The explicit higher twist corrections were 

calculated in Ref. liSlij as expansions in ^ and Fig. 

D contains relative difference between the two approaches 
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and can be seen to be in about permille agreement for 
perturbative values of the energy. 

This remarkable level of agreement implies that in 
order to incorporate one loop massive QCD-flavor thresh- 
old effects into observables which are known only in mass- 
less QCD in the MS scheme, one just has to apply the 
analogous steps as above, namely replacing 

n /- > »AHff(^)= n /.v(^-) (") 

where Q* is the BLM-scale of that process. In order to 
include all flavors one simply has the sum over all m q 
contributions as indicated in Fig. |l|. 

2 Two Loop Results 

At the two loop level the situation becomes much more 
cumbersome. The massless case including pure gluonic 
corrections was calculated in Ref. O. Recently, the mas- 
sive two loop .corrections, depicted in Fig. ||, were ob- 
tained in Ref. 113. 

The double lines denote the static color sources 
for which one employs the heavy quark effective Feyn- 
man rules (HQET), see for example Refs. Ila, while the 
rest contains the full QCD dynamics including massive 
fermion lines. The results in Ref. O were obtained in the 
MS renormalization scheme, related to the MS scheme 
through a simple scale shift 

^ MS = ' (12) 

by a combined analytical and numerical approach. For 
the two point functions a tensor decomposition mas per- 
formed following the techniques given up Ref cB. The 
algebraic manipulation language FORMEJ was used for 
this purpose. The resulting scalar integrals are then 
solved by calculating the pole terms analytically and in- 
tegrating over the remaining Feynman parameters nu- 
merically. 

For the higher point functions this approach is no 
longer applicable due to the presence of the heavy quark 
propagators. It was found to be advantageous to inte- 
grate out the fermion loop analytically and then proceed 
gth the remaining integrations. Details are given in Ref. 

The required expansions in powers of e = 4 — n were 
performed by MAPLE as was the translation of finite 
parts into FORTRAN. This latter point is crucial in face 
of the enormous complexity of the obtained results. 

Strong consistency checks were performed including 
the explicit calculation of all MS counterterms, shown 
in Fig. U to ensure the locality of renormalization con- 
stants, the correct gluon wave function renormalization 
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Figure 3: The two loop massive fermionic corrections to the heavy 
quark potential. The first two rows contain diagrams with a typical 
non-Abelian topology. Double lines denote the heavy quarks, single 
lines the "light" quarks. Color and Lorentz indices are suppressed 
in the first graph. The notation for the remaining digrams is anal- 
ogous. The middle line includes the infra-red divergent "Abelian" 
Feynman diagrams. While the topology of these three diagrams 
is the same as in QED, they contribute to the potential only in 
the non-Abelian theory due to color factors CfCa- In addition, 
although each diagram is infra-red divergent, their sum is infra- 
red finite. The infra-red finite Feynman diagrams with an Abelian 
topology are shown in the last two rows plus diagrams consisting 
of one loop insertions with non-Abelian terms. 



3 



gse,-ct g : gsei-ci f : gse 2 -ct: 




Figure 4: The counterterms corresponding to the massive two loop 
diagrams of Fig. ^. 

constant and agreement with the - pole terms from 
the massless calculation for the remaining three and four 
point two loop functions. 

Furthermore, the absence of infrared divergences in the 
sum of the MS renormalized diagrams M c i, Mvc 3 and 
M-oivc, which all contain infrared poles individually, is 
demonstrated in Fig. |[ denotes the heavy quark 
propagator and only terms with a dependence on fco 
needed to be regulated with a gluon mass. The results 
clearly demonstrate the infrared finiteness of the sum 
which contributes to the physical potential in Eq. [l]. 

The reason why these three amplitudes with an 
Abelian topology do contribute to the QCD potential 
but are absent in QED is connected to the exponentia- 
tion that is implicit in Eq. [y. The terms proportional to 
Cp are actually already included in the potential by the 
exponentiation of the lower order Born and vacuum po- 
larization contributions. In the non- Abelian theory, how- 
ever, there is also a contribution proportional to CfCa, 
which cannot be obtained by the lower order terms and 
thus must be taken into account for the static QCD po- 
tential. 

We also found perfect agreement for the numerical 
results for the finite, MS renormalized expressions in the 
limit m q — > 0. The finite expressions were calculated 
with the Monte Carlo integrator VEGAS Eil and Fig. || 
displays the weighted sum of all two loop MS renormal- 
ized massive fermionic corrections to ay- 
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Figure 5: The sum of the A 2 -dependent amplitudes and countert- 
erms M k c ° + Mv% + M i vc + M k J ct + Mv° 3ct . Circles corre- 
spond to a choice of q 2 = — WGeV 2 and m = m c , triangles to 
q 2 = — lOOGeV 2 and m = m c while the lower curve (squares) has 
q 2 = -lOOGeV 2 and m = m^. the overall normalization neglects 
color factors and the coupling strength. All data are obtained by 
using 10 6 evaluations per iteration with VEGAS and 100 itera- 
tions. The statistical error is indicated and smaller than the sym- 
bols where invisible. The sum for each of the displayed sets of 
parameters is clearly independent of the IR-gluon mass regulator 
A as expected. 
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Figure 6: The complete two loop mass dependence of oty = 

/2 loop 

for m 2 = m 2 = (1.5GeV) 2 (triangles) and m 2 = m 2 = 

(4.5GeV) 2 (open circles). The massless case is also included (line). 
In all three curves we use fi = 0.031. The deviation from the mass- 
less case at the flavor thresholds is of order of 33% and is dominated 
by the new non-Abelian contributions 
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It was found that the overall curve is dominated by 
the non-Abelian threshold behavior (partially due to the 
extra factor of Ca)- The "m c -graph" (triangles) matches 
the massless case for lower values of —q 2 as m 2 <c m|. 
At the respective thresholds we find roughly a 33 % de- 
viation relative to the massless case. This could be very 
significant for applications where quark masses are ex- 
pected to play an important part. At high values of — q 2 
the theory becomes massless and reproduces the leading 
logarithmic terms obtained by the /3-function analysis as 
these coefficients are scheme independent through two 
loops in a massless theory. These analyses can also be 
helpful for the incorporation of massive fermions in lat- 
tice analyses as the heavy quark potential is defined by 
the gauge invariant vacuum expectation value of the Wil- 
son loop in Eq. [l]. 

3 Conclusions 

There is a very elegant and simple way to describe ana- 
lytically massive one loop QCD flavor threshold effects in 
observables calculated in massless QCD in the MS renor- 
malization scheme. All one needs to do is to substitute 
the discontinuous function of active flavors, u/, according 
to Eq. [□]. There is thus no need for complicated higher 
twist calculations in the MS scheme. At two loops, the 
level of perturbation theory to which many observables 
are known in massless (!) QCD, it should be possible to 
extend the one loop analysis presented in section [l] by us- 
ing the now available explicit two loop results discussed 
in section ^. 

The situation is now more complicated, however. 
The only feasible approach is a numerical differentiation 
of the obtained Monte Carlo results. This implies ques- 
tions relating to the technical precision domain. In ad- 
dition one needs to include running mass effects at the 
one loop level as they enter into the-two loop analysis. 
Work towards this end is in progressed and will hopefully 
soon lead to a two loop extension of an effective analytic 
function of quark flavors. 
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